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Abstract. External pressure driven phase transitions of FeSe are predicted 
using ab initio calculations. The calculations reveal that o-FeSe takes transitions 
to NiAs-type, MnP-type, and CsCl-type FeSe. Transitions from NiAs-type to 
MnP-type and CsCl-type FeSe is also predicted. MnP-type FeSe is also found 
to be able to transform to CsCl-type FeSe, which is easier from a-FeSe than 
the transition to MnP-type FeSe, but comparable to the transition from NiAs- 
type FeSe. The calculated electronic structures show that all phases of FeSe are 
metallic, but the ionic interaction between Fe-Se bonds becomes stronger and 
the covalent interaction becomes weaker when the structural phase transition 
occurs from cr-FeSe to the other phases of FeSe. The experimentally observed 
decrease in T c of superconducting a-FeSe at high pressure may be due to a 
structural/magnetic instability, which exists at high pressure. The results suggest 
us to increase the T c of ct-FeSe if such phase transitions are frustrated by suitable 
methods. 
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1. Introduction 

FeSe has been studied for its magnetic properties [U El H] and showed a phase 
transition from a-FeSe to /3-FeSe (tetragonal to hexagonal) as the growth temperatures 
were increased [2J. Interest in a-FeSe have grown after Hu et al. [I] observed 
superconductivity with a T c ~ 8 K because it is structurally related to the newly 
existing investigations of FeAs superconductors [5] and subsequently the T c was 
increased by the application of pressure [6, 7\ and doping [8] . External pressure up to 
2.2 GPa enhances both spin fluctuations and T c in a-FeSe [5]; however a large external 
pressure further decreases T c [5J. 

Experimental results show that tetragonal FeSe undergoes a structural transition 
to a hexagonal phase [TU] and at 12.5 GPa to an orthorhombic phase [IT]. More 
recently, Medvedev et al. [& reported the increase of T c from 8.5 to 36.7 K under 
an applied pressure of 8.9 GPa and also showed a phase transition from tetragonal 
to NiAs-type hexagonal at 12 GPa. Compared with the FeAs superconductors [5], 
which also show a structural phase transition from tetragonal to orthorhombic [T2lll3j . 
a-FeSe has not only the same planar sublattices, but also displays structural and 
magnetic instabilities [14]. 

Albeit FeSe shows structural transitions, it is not clear how the electronic- 
structure responds to external pressure, and what the transition path of FeSe should 
be. It is then invaluable to study FeSe in different crystallographic environments 
to observe the effects of structural changes upon changing the atomic volume or 
applying pressure. Thus, because of its comparative chemical simplicity with FeAs 
superconductors, FeSe provides a unique opportunity to study the interplay of the 
structure, magnetism, and superconductivity. Therefore, the determination of the 
structural and magnetic phase transitions in FeSe is essential for understanding their 
electronic properties. To our knowledge, no systematic study has been performed 
which can explore the mechanism of the structural phase transitions in stoichiometric 
FeSe. We show that there are structural and magnetic phase transitions at high 
pressure. Such phase transitions also affect the magnetic, electronic, and bonding 
properties of FeSe; these are considered as an indirect clue for improving the T c of 
FeSe. 

This paper is organized as follows. We describe the computational models and 
methods in Sec.[2]which are used to investigate the structural phase transitions of FeSe. 
In Sec. [3] we discuss our calculated results, i.e., structural phase transition, magnetism, 
and electronic structures of FeSe. The results and discussions are summarized in Sec. [4] 

2. Computation Models and Methods 

To search for the structural phase transition of FeSe, we considered different types 
of crystal structures, i.e., anti-PbO — commonly known as a-FeSe, as well as NiAs-, 
MnP-, CsCl-, and CuAu-type FeSe (see Fig. [lj. These structures are related with 
each other, e.g., MnP-type FeSe can be considered as a deformed NiAs-type FeSe, 
changing the local hexagonal symmetry around the Fe atoms in the NiAs-type phase 
to tetrahedral results in the a-FeSe phase. The detailed structure relationship will be 
discussed in Sec. |3] 

Ab initio calculations were performed using the total-energy all-electron full- 
potential linearized augmented plane-wave (FLAPW) method [TS] in the QMD- 
FLAPW package [TB] based on both the local spin density approximation (LSDA) [T7] 
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and the generalized gradient approximation (GGA) [18]. Integrations inside the 
Brillouin zone (BZ) were performed by the improved tetrahedron method [TS] over 
a 15 x 15 x 15 mesh within the three dimensional (3D) BZ and with an energy cutoff 
at 4.10(27r/et), where a is the lattice parameter of each calculation. A 16.36 (2ir/a) 
star-function cutoff was used for depicting the charge and potential in the interstitial 
regions. Lattice harmonics with I < 8 were employed to expand the charge density, 
potential, and wave-functions inside each muffin-tin (MT) sphere of radii 2.2 a.u. for 
Fe and 1.9 a.u. for Se. The convergence of the computation parameters were checked 
carefully [20] • 

All core electrons were treated fully relativistically and valence states were 
calculated scalar relativistically, i.e., without spin-orbit coupling 21J. The explicit 
orthogonalization (XO) procedure was employed to ensure the orthogonality between 
the semicore and valence states with the spill-out of each MT sphere semicore charges 
being superposed [22 . Self-consistency was assumed when the difference between 
input and output charge densities became less than 1.0 x 10 -4 electrons/a. u. 3 

The calculated total energy-volume (E-V) data were fitted to the Birch- 
Murnaghan equation of states 23J which enables us to calculate the enthalpy H = 
E + pV, where E is the internal energy, p is the external pressure, and V is the 
volume of the system. All calculations were carried out for nonmagnetic (NM) and 
ferromagnetic (FM) states at different volumes. The optimized lattice parameters of a- 
FeSe were calculated to be a = 3.747 (3.701) A, c/a = 1.474 (1.505), and volume/atom 
ratio = 19.685 (19.078) A 3 with the GGA (LSDA) results which are comparable to the 
recent experimental values jllj . The calculated internal coordinates and the formula 
unit volumes of the considered crystals at individual equilibria are summarized in 
Table [T] We carefully checked that the physics of these systems does not depend on 
the details of internal parameter values. 

3. Results and Discussions 

3.1. Structural Phase Transition 

In Fig. [2] one can see that a-FeSe is the most stable structure and under sufficient 
volume contraction (external pressure), a-FeSe becomes unfavorable compared to the 
other phases of FeSe. Around the equilibrium volume of a-FeSe all other phases are 
in metastable states by an energy barrier per formula unit (f.u.) of order 1.0 eV. As 
we compress the unit cell volume, a-FeSe goes to a higher energy state and NiAs-type 
FeSe becomes a stablized phase. NM a-FeSe starts to cross FM NiAs-type FeSe at 
a volume of ~ 30.74 A 3 . Further compression transits a- FeSe to the NM NiAs-type 
phase at 30.19 A 3 , but before crossing the NM NiAs-type FeSe phase, a-FeSe crosses 
the energy curve of MnP-type FeSe. Hence, a slight compression deforms NiAs-type 
FeSe to MnP-type FeSe, where the hexagonal symmetry is broken and can be viewed 
as a distorted NiAs structure. We can also see that NM a-FeSe begins to overcome 
FM MnP-type FeSe at about 30.33 A 3 and NM MnP-type FeSe at about 29.86 A 3 . 

NiAs-type FeSe can be transformed to CsCl-type FeSe if we push the Se atom 
to the center of the unit cell and simultaneously deform the a and b axes of the 
hexagonal lattice to a cubic lattice under a large external pressure, to break the bonds 
between Fe and Se atoms in the hexagonal structure for making new bonds in the 
cubic structure. In other words, tetragonal a-FeSe can be thought of as a distorted 
CsCl-type FeSe with an elongated c axis. Hence, one can naturally expect a transition 
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to CsCl-type FeSe at high pressure corresponding. On the other hand, MnP-type 
FeSe is a distorted structure where phase transition to CsCl-type FeSe will require 
larger pressure as compare to NiAs-type FeSe. Although CuAu-type FeSe lies much 
higher in energy, extrapolating a-FeSe to a small volume region tells us that a-FeSe 
can enter into a metastable region at 27.80 A 3 and CuAu-type FeSe may become a 
stable structure. However, the common tangent does not pass through both curves 
which makes difficult the a-FeSe to CuAu-type FeSe transition. We must note that 
Fig. [2] can give an indication of structural phase transitions of a- FeSe under pressure, 
and the volumes discussed above are the critical volumes at which a-FeSe crosses the 
binding curve of the other phases, e.g. NiAs-type FeSe 

The above discussion revealed the relative structural stability of a-FeSe. However, 
to have complete understanding of the phase transition of a-FeSe, we calculated the 
transition pressures. The transition pressures (p t ) from a- FeSe to the other structure, 
e.g., NiAs-type FeSe, can be obtained either from the slope of the common tangent 
of both E-V curves in Fig. [2] or from the usual condition of equal enthalpies, i.e., 
the pressure pt at which the enthalpies H of a-FeSe and NiAs-type FeSe are the 
same. We follow the latter approach, since it is well known that the thermodynamic 
stable phase at some given pressure and zero temperature is that one with the 
minimum enthalpy |24j . At p t , two phases have the same enthalpy, and then the 
transition pressure is determined by equating the enthalpies of the two phases of 
FeSe. Note that this is a standard procedure to calculated the transition pressure of 
materials [25j [26l [27] . To see the possible phase transitions, Fig.|3]shows the calculated 
enthalpy curve of a-FeSe, NiAs-, and MnP-type FeSe as a function of p. 

Following the above procedure, the transition pressures of FeSe are calculated and 
the results are summarized in Table [2j The GGA (LSDA) p t of a- FeSe to NiAs-type 
FeSe is calculated to be ~ 17.5 (13.0) GPa, whereas pt of a-FeSe to MnP-type FeSe is 
~ 17.70 (13.05) GPa which are in agreement with recent experimental reports [6l 111], 
It is clear that a slight external pressure deforms the hexagonal structure to the 
orthorhombic one. The experimental reports state that structural phase transition 
from a-FeSe to MnP-type requires ~ 0.5 GPa more extra pressure than does NiAs- 
type [SI [TT]. Our calculations follow the same trend and the GGA (LSDA) results 
demonstrate that ~ 0.20 (0.05) GPa extra pressure than NiAs-type is required for 
the MnP-type FeSe. We also calulated the transition pressure of a-FeSe to CsCl- 
type FeSe which is ~ 20.45 GPa. Therefore, our calculations reveal that a-FeSe takes 
transition to NiAs-type FeSe, MnP-type FeSe, and CsCl-type FeSe at high pressure 
irrespective of their magnetic structures. The predictions are in excellent agreement 
with recent experiments O [TT] . Table [2] also shows that there is a possibility of a 
phase transition from NiAs-type to CsCl-type FeSe at ~ 42 GPa. We can see that it 
is easier for a-FeSe to take transition to NiAs- and then to MnP-type FeSe. However, 
large external pressure ~ 60 GPa is required for MnP-type to CsCl-type FeSe phase 
transition. Further recent experimental observation found CsCl-type FeS [28 , which is 
isoelectronic to FeSe, at high pressure and this supports the expectation of CsCl-type 
FeSe [29]. 

3.2. Magnetism 

The calculated magnetic moments per Fe atom (//pe) are shown in Fig. [4] We did not 
find any magnetism in a-FeSe in the whole volume range, which is in agreement with 
recent first-principles calculations [3U] and experiment [2] . The steps in fip e of NiAs- 
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type FeSe indicate that the magnetism of FeSe is sensitive to volume compression and 
expansion. Below the volume of 22 A 3 , NiAs-type FeSe has zero magnetic moment 
and the onset of magnetization can be seen for volumes greater than 22 A 3 ; the small 
step in yitFe is due to the onset of magnetization. Once the magnetization develops, 
/iFc increases slightly and has a constant slope in the volume region 23-25 A 3 ; this 
is the region where MnP-type FeSe is more stable in energy than NiAs-type FeSe. 
The step in /ip e happens at ~ 31 A 3 , because near this volume the hexagonal FeSe 
is more stable in energy than tetragonal FeSe. Interestingly similar steps in /xp e can 
also be seen for MnP-type FeSe, where the onset of magnetization shifts to a higher 
volume region as compared with NiAs-type FeSe, i.e., the onset of magnetization 
starts at ~ 25 A 3 , exactly at this volume that the FM MnP-type enters into FM 
NiAs-type FeSe. Then /ip c increases with increasing lattice volume and the change in 
slope can be seen around the volume region where a-FeSe crosses the energy curve 
of MnP-type FeSe. The CuAu-type FeSe has zero //p e for volumes < 22 A 3 and /iF e 
monotonically increases with the lattice volume. CsCl-type FeSe is more interesting 
due to its nonzero fj,p e i n the whole region of volume and this may indicate that while 
even all the phases of FeSe become NM, CsCl-type FeSe retains its magnetic structure 
in an unconventional manner [29] . 

3.3. Electronic Structures 

The phase transitions mentioned are understood by comparing the electronic 
structures in terms of the density of states (DOS). The calculated total and atom 
projected DOS in the FM states are shown in Fig. [5| The total DOS of a- FeSe shows 
metallic behavior which is in agreement with recent photoemission and transport 
measurements [I]. The exchange splitting is absent in a-FeSe and the Se p-states 
lie well below the Fermi energy (Ep), The Se p-states are hybridized with the Fe 
d-states and this hybridization is strong at a lower energy, ~ —4.0 eV. The strong 
peak near E-p is formed mainly from the Fe e s -states, and the Fe ^g-states are more 
delocalized in energy and extend over a wide range of energy. The DOS shows that 
Ep is dominated by the Fe d-states, which is also found in recent experimental results 
on a-FeSe [3T] and the FeAs superconductor [35]. One can also see a pseudogap in 
the electronic structure of a-FeSe which was also found in previous calculations [33 
and in the FeAs based superconductor [34]. This reveals that the electronic structure 
of a-FeSe is essentially similar to those of FeAs based superconductors. 

The NiAs-type FeSe DOS in Fig. [5](b) shows a metallic property and now there is 
an exchange splitting between the majority and minority spins, which was absent in 
a-FeSe. Therefore, the structural phase transition is accompanied by a magnetic phase 
transition. Simultaneously the Fe local symmetry is also changed from tetrahedral to 
octahedral. The Se p-states extend to lower energies and strongly hybridize mainly 
with the Fe i2g-states. Furthermore, now there is no pseudogap in the majority spins, 
but it exists in the minority spin states. The peak at ~ —2.0 eV is contributed by the 
Fe i2g-states. This switch of e g to ti g states is due to the change of local symmetry. 
Once the local symmetry of the Fe atom is changed to hexagonal (octahedral) , further 
breaking of hexagonal symmetry will not have much affect on the electronic structure 
as seen in MnP-type FeSe DOS in Fig. pjfc). The electronic structure is identical to 
NiAs-type FeSe except for some minor changes around Ep in the majority and minority 
spins. In the majority spins, Ep lies in the valley; in the case of minority spins, the 
splitting between the t2 g and e g states at Ep is increased and this is probably due to 
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the breaking of hexagonal symmetry. 

In the case of NiAs-type FeSe, these minority spin states were nearly degenerate, 
but the lattice distortion removes this degeneracy and the e s -states move to higher 
energy in MnP-type FeSe, lower the electronic energy and E-p falls in the valley. One 
can see that the lattice distortion significantly decreased the DOS at E-p, and the 
MnP-type FeSe has a lower DOS at Ep than does NiAs-type FeSe. In the case of 
CsCl-type FeSe, the majority Fe e g -states are completely occupied and the ^-states 
are partially occupied. Metallic conduction is mainly taken by the ti g spins in the 
majority as well as in the minority spin states. The total DOS shows the mixing 
of the Fe-d and Se-p states. The total DOS of CuAu-type FeSe also shows metallic 
character; the hybridization between the Fe-d and Se-p electrons can be seen. Note 
that both the t^g and e g spins are completely occupied in the majority states. The 
exchange splitting between the t 2g and e g spins in the majority state is weaker than 
in the minority spin states. As we saw that our calculated systems are metallic in in 
low and high pressure regions and the electronic structures agree with the previous 
experimental and theoretical data which suggest us that the effect of onsite Coulomb 
interaction, i.e., LDA+U type interaction in FeSe will not have significant affects. 

To shed more light on the structural phase transition of a-FeSe and to see how 
the chemical bonding changes, the charge densities were also calculated and plotted 
in Fig. [6] in a plane that contains both the Fe and Se atoms. The a-FeSe phase forms 
an enhanced density on the side of the Fe atom. It is clear that the a-FeSe bonds 
are more covalent, a bonding that causes the apparent deviation of the charge density 
contours from a spherical shape. Inside the Fe-Se block, the metal-like Fe-Fe (Se-Se) 
bonding occurs due to delocalized Fe 3d (Se p) states. Interestingly, similar bonding 
was also found in the FeAs superconductors [35]- Note that such a bonding is either 
absent or weak in the other crystal structures of FeSe. The structural transition to 
NiAs-type FeSe also changes the bonding character between the Fe and Se atoms. 
The charge density of NiAs-type FeSe, see Fig. [6](b) , clearly shows that the bonding 
between the Fe and Se atoms is less covalent and the charges are more localized at 
Fe and Se atoms and are nearly spherical. The bonding in MnP-type FeSe is similar 
to the NiAs-type FeSe. In the CsCl-type FeSe and CuAu-type FeSe, the bonding 
becomes more ionic and the charge density becomes nearly spherical. It is obvious 
that the ionic interaction between the Fe-Se bonds becomes stronger and the covalent 
interaction becomes weaker, when the structural phase transition occurs from a-FeSe 
to the other phases of FeSe. We must note that the role of defects (Se or Fe vacancies) 
in FeSe can not be ignored and such defects, which are usually in high energy states, 
may change the magnetic structure of FeSe [3D]. The effect of Se/Fe vacancies on the 
structural phase transition of FeSe is beyond the scope of the present work. 

There might be a relation between superconductivity and covalency because a- 
FeSe is more covalent than the other phases of FeSe, which have not yet shown 
superconductivity. The ground state structure of FeSe is NM a-FeSe and it does 
not show any magnetic or structural instability at zero pressure. The experimentally 
observed decrease in the T c of a-FeSe at high pressure [6] might be due to a 
structural/magnetic instability because our calculations show that there is a structural 
and magnetic instability under high pressure and a structural phase transition takes 
place to non-superconducting FeSe, e.g., NiAs-type FeSe. Although their is no 
concensus on the mechanism of the observed superconductivity of a-FeSe, the Fermi 
surface topology plays a crucial role [35] for forming the mediating quasi-particles for 
Cooper pairs as in cuprates [37]. Once a-FeSe transformed into the other structures, 
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such Fermi surface topology and the corresponding van Hove singlularity near the 
Fermi surface will not be maintained anymore. Therefore, the decrease in T c at high 
pressure may be due a structural/magnetic instability. This suggests that we may 
increase the T c of FeSe, if such a structural phase transition at high pressure is avoided. 
It may be suggested that the structural phase transition in a-FeSe can be controlled by 
doping [38, 39 . Further work is necessary to control this phase transition or increase 
the transition pressure of a-FeSe. 

4. Summary 

Using ab initio calculations, we showed the phase transition of superconducting a- 
FeSe. The predicted structural phase transition and transition pressure were shown 
to be in excellent agreement with existing experimental observations. Our calcula- 
tions also showed that there is also a structural phase transition from NiAs-type to 
MnP-type and CsCl-type FeSe. Phase transition from MnP-type to CsCl-type FeSe 
requires larger pressure than NiAs-type FeSe. The electronic structure and charge 
density were similar to the FeAs based superconductors. The calculations revealed 
that the covalent bond between Fe and Se becomes weaker and the ionic bond be- 
come stronger when a-FeSe is transformed to other structures. The experimentally 
observed decrease in T c at high pressure of a-FeSe may be due to a structure/magnetic 
instability. A possible scenario for the enhancement of T c was also discussed. 
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Table 1. Structural parameters of FeSe. The last column shows the optimized 
volume V (in units of A 3 ) per formula unit 



Type 
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Atomic positions 
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a— FeSe 


Fe(2a) 


; ± ±0 
\ u ; g u 
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Se(2c) 


0.243 






NiAs-type 


Fe(2a) 
Se(2b) 


0; ± 

12 21 ,1 

3 3 " ' 3 ' 3 ' " ' 2 


0.250 




33.640 


CuAu-type 


Fe(2a) 
Se(2b) 


; ±±0 
ini - nn 

2 2 ' 2 2 






31.890 


MnP-type 


Fe(4c) 
Se(4c) 


±(" |; 5 ~U,V+ |, j) 

±(" J! 2 + 2' 4) 


0.200 
0.570 


0.005 
0.190 


29.000 


CsCl-type 


Fe 
Se 



111 

2 2 2 






28.370 



Table 2. Calculate structural transition pressures (in units of GPa) of FeSe. The 
available experimental data are taken from Ref. |11| [6). 
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Experiment 
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MnP-type 
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CsCl-type 
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MnP-type-^ 


CsCl-type 
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Tables 
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Figure 1. (Color online) Crystal structures of FcSc.(a) a- FeSe, (b) NiAs-, (c) 
MnP-, (d) CsCl-, and (e) CuAu-type FeSe. Red and yellow balls show the Fe and 
Se atoms, respectively. 



Figure 2. (Color online) The calculated total energy in units of eV/f.u. as a 
function of lattice volume (A 3 ) for ct-FeSc (a) NiAs-, (b) MnP-, (c) CsCl-, and (d) 
CuAu-type FeSe. The total energy was measured with respect to the total energy 
of a-FeSe equilibrium volume. Filled circles (squares) show the total energy in 
the NM (FM) state. Red, black, cyan, green, and blue colors represent ct-FcSe, 
NiAs-, MnP-, CsCl-, and CuAu-type FeSe, respectively. The inset (e) shows the 
E — V curve of all the structure in volume interval 20-35 A 3 . 



Figure 3. (Color online) The calculated enthalpy in units of Ry as a function of 
pressure (GPa) for (a) a-FeSe — > NiAs-typc, (b) a-FeSe — > MnP-type, (c) a-FeSe 
-> CsCl-typc, (d) NiAs-type -s> CsCl-type, (e)NiAs-type-s> MnP-type, and (f) 
MnP-type — > CsCl-type FeSe. Red, black, cyan, and green colors represent a- 
FeSe, NiAs-, MnP-FeSe, and CsCl-type FeSe respectively. Dotted lines represent 
the NM states 



Figure 4. (Color online) Calculated magnetic moments (in units of fi-g) per Fe 
atom as a function of lattice volume (in units of A 3 ) of FM FeSe. Red, black, 
cyan, green, and blue colors represent a-FeSe, NiAs-, MnP-, CsCl-, and CuAu- 
type FeSe, respectively. 

Figure 5. (Color online) Total and atom-projected DOS in the FM states of 
(a) a-FcSc, (b) NiAs-, (c) MnP-, (d) CsCl-, and (c) CuAu-type FeSe at their 
equilibrium lattice volumes. Red, blue, and green lines represent Fe t2 g , Fe e g 
and Se p states, respectively, whereas black lines in the bottom panels represent 
the total DOS. The Se p states are multiplied by a factor of 10 and the Fermi 
levels (E-p) are set to zero. The total DOS is given per formula unit for comparison 
purpose. 

Figure 6. (Color online) Charge-density contours of (a) a-FeSe,(b) NiAs-, (c) 
MnP-, (d) CsCl-, and (e) CuAu-type FeSe in the planes containing Fe and Se 
atoms. The lowest contour starts from 2x 10 -4 electrons/a.u. 3 and the subsequent 
lines differ by a factor y/2. The color code of the charge density contours is ordered 
as green, yellow, and red, in ascending order. The dotted lines show the bond 
between Fe and Sc atoms. 



Figure captions 



Fig. 1(a) 




Fig. 1(c) 



Fig. 1(d) 




Fig. 1(e) 




I I I I I I 

15 20 25 30 35 40 45 50 
Volume (A 3 )/f.u. 



Fig. 2(a) 



NM-MnP 
FM-MnP 




i i i i i i 

15 20 25 30 35 40 45 50 
Volume (A 3 )/f.u. 



Fig. 2(b) 




I I I I I I 

15 20 25 30 35 40 45 50 
Volume (A 3 )/f.u. 



Fig. 2(c) 



19 



3 

2.5 ^ 

2 
1.5 
1 

0.5 





NM-CuAu • 

FM-CuAu ■ 

NM-PbO • 
PbO 



15 20 25 30 35 40 45 50 
Volume (A 3 )/f.u. 



Fig. 2(d) 




Fig. 2(e) 




-0.2 1 1 1 1 1 1 1 

5 10 15 20 25 30 35 4 

Pressure (GPa) 



Fig. 3(a) 




Fig. 3(b) 




-0.2 1 1 1 1 1 1 1 

5 10 15 20 25 30 35 40 

Pressure (GPa) 



Fig. 3(c) 



24 




Pressure (GPa) 



Fig. 3(d) 



25 




Pressure (GPa) 



Fig. 3(e) 



20 




0.4 1 - — — — — — — -J 

40 45 50 55 60 65 70 75 80 



Pressure (GPa) 



Fig. 3(f) 



27 



4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 


-0.5 



FM-MnP 
FM-NiAs 
FM-CsCl 
FM-NiAs 
FM-CuAu 




: 

■ 


- 


y \ 


■ 

■ 

■ 

■ 

■ 

■ 


f ■ 

n 





15 



20 



25 



30 



35 



40 



45 



50 



Volume (A°) 3 /f.u. 



Fig. 4 



28 




Fig. 5(a) 



29 




Fig. 5(b) 



30 




Fig. 5(c) 




Fig. 5(d) 




Fig. 5(e) 




Fig. 6(a) 




Fig. 6(b) 




Fig. 6(c) 



Fig. 6(d) 




Fig. 6(e) 



